Electromagnetic turbulence driven by the mixed mode instability 
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In continuation of previous work, numerical results are presented, concerning relativistically 
counter-streaming plasmas. Here, the relativistic mixed mode instability evolves through, and be- 
yond, the linear saturation - well into the nonlinear regime. Besides confirming earlier findings, that 
wave power initially peaks on the mixed mode branch, it is observed that, during late time evolution 
wave power is transferred to other wave numbers. It is argued that the isotropization of power in 
wavenumber space may be a consequence of weak turbulence. Further, some modifications to the 
ideal weak turbulence limit is observed. Development of almost isotropic predominantly electro- 
static - partially electromagnetic - turbulent spectra holds relevance when considering the spectral 
emission signatures of the plasma, namely bremsstrahlung, respectively magneto-bremsstrahlung 
(synchrotron radiation and jitter radiation) from relativistic shocks in astrophysical jets and shocks 
from gamma-ray bursts and active galactic nuclei. 



Counter-streaming plasmas subject to the relativis- 
tic mixed mode instability (MMI), driven by relativistic 
beams of electrons and ions, have previously been exam- 
ined theoretically [j], Q , and numerically with particle-in- 
cell (PIC) simulations (3,0]. Growth rates of the mixed 
mode, the filamentation and the two-stream instabilities 
(hereafter MMI, FI and TSI, resp.), are 

7mm = \/3(n 6 /16n ?) ) 1 / 3 r(i; b )- 1 / 3 

7 f = (Vc)(n6K) V2 I>6)~ 1/2 
7 T5 = V3(n b /16n ?) ) 1 / 3 r(i; b )- 1 , 



with rib, n p , Vb and T(vb) = (1 



fj/c 2 ) the beam 



and background plasma density, the beam speed and 
beam bulk Lorentz factor, respectively, in the back- 
ground rest frame. For a large volume in parameter 
space, V € {rib, n p, v b, T(vb)}, the MMI will have the 
highest growth rate of the three possible instabilities for 
the relativistic beam-plasma interaction. The MMI prop- 
agates at oblique angles with respect to the beam velocity 
^(kMMJjVb) oc arctan [(wb/wp) 1 / 2 ], where v p is the back- 
ground plasma thermal speed and v p <C Vb A thor- 
ough theoretical discussion of the relationship between 
the MMI, FI and TSI modes is given in 0. 

Due to its mixed nature, the MMI contains both an 
electrostatic and an electromagnetic wave component Q . 
Potentially, both electrostatic and electromagnetic tur- 
bulence (wave mode coupling leading to cascades/inverse 
cascades in fc-space) is possible in such systems. This po- 
tential for producing very broad band plasma turbulence 
(in both E and B fields) is highly relevant when consider- 
ingfor example inertial confinement fusion experiments 
[ETUI . Other examples where electromagntic wave turbu- 
lence is important are astrophysical jets and shocks from 
gamma-ray bursts and active galactic nuclei, where ambi- 
ent plasma streams through a shock interface moving at 



relativistic speeds - see e.g. @, H, El- Further studying of 
the MMI in particular, here, beyond the linear regime to 
the development of turbulence is thus highly motivated. 

Performing high resolution 2.5D PIC code simulations, 
we examine the turbulent wave spectra, that develop dur- 
ingsaturation of the MMI by the trapping of electrons 
(JHI- The code solves Lorentz' equations of motion for 
an ensemble of computational macro-particles (CPs) de- 
fined in 2D3V; {x,y,p x ,p y ,p z ,t}, and Maxwell's equa- 
tions in two spatial dimensions; {Ei x . y , B 2 }(i£, y, t). A 
PIC code scaling is chosen wherein m e = e = c = 1 for 
the electron mass, elementary charge and light speed. 

We model a relativistic beam consisting of co-moving 
electrons and protons with a mass ratio mi/m e = 1836 
and equal number densities. The background plasma 
consists also of equal number densities of electrons and 
protons. The computational system is consequently 
charge neutral and current neutral and neither elec- 
tric, nor magnetic, fields are present, initially. The 
background plasma is chosen as our laboratory frame. 
Here, the bulk flow velocity vector of beam electrons 
and beam protons is set to v b = u b y w 0.97c y, i.e. 
with T(vb) = 4. Densities and temperatures, for the 
cold beam and warm background plasmas are set to 



rib 



0.1 and T b = 10" 4 (v th ,b = 0.01c), and n v 



0.9 



and T p = 10 



2 {vth,p = 0.1c), using the aforementioned 
scaling of the relevant natural constants. With these 
parameter choices the plasma frequency is Lu pe ps 0.9, 
and the electron skin depth becomes S e = c/uj pe 1.1. 
Time is henceforth normalized to the plasma frequency, 
t = t'uj pei and wavenumbers to the volume dimensions, 
i.e. k x = 1 = 2ir/L x . All boundary conditions on both 
the fields and particles are periodic. 

To ensure robustness and numerical accuracy of our 
results we have compared two simulations with the same 
physical parameters but with different grid and parti- 
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FIG. 1: (Color) (a) The electric and magnetic energy densities eg (red upper curves) and es (blue lower curves) measured for 
Simulation 1 (dashed curves) and Simulation 2 (solid curves). Noise builds up in the region 1 and the MMI develops linearly 
in the region 2. The ratio e^/es increases after the saturation, at iw e ,i ~ 50. We evaluate the time evolution of wave power 
spectra on intervals A through D. (b) the power spectrum, P(k x , k v ,t = 36), corresponds to spectrum of the MMI during linear 
growth, (c) shows the same spectrum just after saturation, P(k x ,k y ,t = 52). k x , y are in units of 2Tv/L x ,y The color scale is 
10-logarithmic and identical for both (b) and (c). 



cle resolutions. Simulation 1 uses a grid size L x = 
1500AZ = 74<5 e , L y = 1500A y = 74<S e where A x = A y 
is the grid spacing. The plasma is resolved in particle 
phase space with N e = N p = 200 CPs/cell/species - a 
total of N t = 400 CPs/cell. In Simulation 2 the cor- 
responding run parameters are: L x = 59<5 e = 1200A,,;, 
L y = 59<5 e = 1200A y , A^ = A y and N t = 600 CPs/cell. 
Figure QJa) displays the electric and magnetic energy 
densities, integrated over the simulation volume. The 
data are normalized to the total electric field energy for 
Simulation 1 at t = t en d = HOo;^, 1 . As seen from fig- 
ure [lja) we obtain excellent agreement for the develop- 
ment of the energy densities through all stages of the 
MMI instability in the reference run. The small discrep- 
ancy in energy is solely due to a slightly higher initial 
noise level of Simulation 2 that persists throughout; this, 
in turn, shows us that energy conservation is also agree- 
able. Having checked that all other results presented here 
are also robustly reproduced in Simulation 2, we present 
data solely from Simulation 1. We avoid re-examination 
of the electron distributions and refer to earlier work 
treating the linear MMI phase [H, 0] and dynamics of the 
ions. We merely note here that the ions begin to undergo 
instability in the late linear MMI phase and has impor- 
tant influence on long time evolution. The remainder of 
this brief communication focuses on the field dynamics 
during transition to non-linear development. 

Figure [3a ) displays the total field energy history. Ini- 
tially, fields are due to noise and practically electrostatic 
already after t ~ 5 (Region 1 in FigJTJa) ). Magnetic 
fluctuations, however, rapidly build up to equipatition 
with the electric fluctuations. Once in equilibrium, the 
ratio Eb grows in unison with Ee through the late linear 
stage (Region 2) of the MMI, albeit at a slightly lower 
value. This ratio again increases after MMI saturation 



at t > 50. It reaches a quasi-steady value Ee/Eb ~ 10 
at the simulation's end (60 < t < t en d = HO), signifying 
the predominance of electrostatic fields. 

Figure DJb,c) displays power spectra PE(k x ,k y ,t) = 
\FFT[E x (t) + iE y (t)}\ 2 of the electric field, in the quad- 
rant {k x ,k y } > {0,0} at times t = 36 (during linear 
MMI) and t — 52 (just prior to saturation). We have 
rendered only a single quadrant for the sake of clar- 
ity and display resolution. The spectral representation 
of the electric field, PE{k x ,k y ,t) reveals the unstable 
wave branch, and the k-range of the turbulent spec- 
trum. With 5 e ~ 1.1 and v b w 0.97c (from T(v b ) = 4), 
we find the most unstable mode at k y w (0.98<5 e ) _1 or 
k y = ky/(2ir /74<5 e ) sa 11, which agrees well the k y of 
the MMI branch in FigQJb,c). Significant wave power is 
present over a broadened range of k x ,k y at tu) e> i = 52, 
which peaks on the initial MMI branch at t < 36. The 
magnetic power spectrum PB(k x ,k y ,t) = \FFT[B z (t)]\ 2 
(not shown) - also due to the MMI [2| - follows closely 
that of the electric field albeit at lower power until about 
t r-j 50 when linear phase has saturated. Also, power is 
slightly enhanced at large k x , k y due to the finite grid 
instability, which was treated in Q . 

Figure [T{b,c) also reveals a rapid time dependent 
spreading of wave power away from the mixed mode 
branch. This effect is not due to Larmor effects of the 
magnetic field; even if assuming a peak magnetic field 
B m ax ~ 0.2 (measured) everywhere during 20 < t < 60, 
yields a relativistic electron gyro frequency of uo ce ~ 
Lo pe (eB max ) I ' (m e Y(vb)c) ~ 0.05cj pe . Hence, the Larmor 
time scale is not a governing time scale. Langmuir wave 
turbulence has been found [10( to spread into a ring dis- 
tribution power spectrum in k-space, assuming a weak 
turbulence model. This result was obtained for initially 
mono-directional Langmuir electrostatic turbulence, in 
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FIG. 2: (Color) Temporal development of the power spectrum ratio, R(k x ,k y ,t, At), at the times tiui e ,i = 36, 44, 52, 
60 - (a), (b), (c) and (d) respectively. Color scaling is 10-logarithmic, and zero-level (blue) equals zero growth rate; i.e. 
Logio(-R(fc^, k y , t, At)) = 0, following from a ratio P{tf) / Pfa) = 1. Only one quarter, {1 < k x < 375, 1 < k y < 375}, of the 
full spectrum quadrant, for {k x > 0, k y > 0}, is shown for purposes of visual clarity. Note the fc-scale is linear. 



the non-relativistic limit [icj |. Our simulations - on the 
other hand - fully resolve all particle and wave dynam- 
ics. Consequently, we capture the spectral dynamics of 
both weak turbulence (wave- wave interaction) and strong 
turbulence (wave-particle interaction). By analysing the 
rapid temporal evolution of PE(k x ,k y ,t), we can deter- 
mine whether the observed spread of the power in k-space 
is partially or fully due to weak Langmuir turbulence. 
The k-interval with the strongest change should rotate 
around k x ,k y = on shells of |fc| = const. The quantity, 
of interest is 

R{k x , ky, h, A t ) = log 10 (P E [k x ,ky,ti]/PE[k x ,k y ,to\) , 

which captures the difference of PE(k x ,k y ,t), at 
two successive simulation time steps. We calculate 
R(k x , k y , ti, A t ) at four times t± = 36, 44, 52 and 60, with 
At — (t\ — to) = 8, corresponding to the time intervals 
"A" through "D" in FigOJa). 

Indeed, figure [2] reveals rich structure and dynamics 
of R(k x , k y , ii, At). Initially (t < 36) power grows uni- 
formly, covering a range, 1 < k x < 280, and 1 < k y < 50, 
mainly along fc^-axis. Later it rotates away from the 
fc^-axis, towards the axis fc y -axis while maintaining its 
elongated shape. Power is transferred between the x- 
and y-directions on curves of approximately constant |fc|. 
We conclude that, qualitatively, what is observed from 
FigEIa-d) is a consequence of driven weak turbulence. 
Some modification to pure (fe £ ,fc y )-rotation results from 
relaxing the assumptions of [lOfl . A much weaker struc- 
ture disperses (both positive (light green) and negative 
(blue) growth) wave power produced on the finite grid in- 
stability 0] branch. It is eventually outpowered by MMI 
driven turbulence. 

Ideal weak turbulence scatters wave energy [l(| on 
shells azimuthally centered in /c-space (k x ,k y ). On the 
other hand strong turbulence, mediated by wave-particle 
interactions and phase space hole coalescence, trans- 
fers wave energy to lower wavenumbers along k y , hence 



perpendicular to the wave scattering direction of weak 
turbulence (which occurs on constant fc-shells). Wave- 
particle interactions in mono-directional turbulence were 
considered previously 11 , [13] , as was the interplay be- 
tween electrons and isotropic Langmuir turbulence [HI]. 
The MMI in our results, however, furthermore yields 
turbulent magnetic fields, which were not taken into 
account by 11, 12, [lj| ■ At late times the magnetic 
field will participate increasingly in producing particle 
trapping and contribute magneto-bremsstrahlung (syn- 
chrotron/jitter Q) to the radiation signature of the 
plasma. 

Lastly, figure [3fa-c) displays azimuthally integrated 
power spectra of the electric field for 10 different 
timesteps during the simulation. Furthermore, for a 
single time, t — t en d = 110, and Fig(3]^d) shows the 
corresponding spectrum for the magnetic field. At the 
earliest times (low power) we observe, again, growth 
of small scale noise which is seen to be approximately 
thermal. This noise is generic to PIC codes. However, 
Figj3jb), P(E x (k ( j > ))), we see that the high k noise 
grows further at all times. This we attribute mainly to 
the finite grid instability in the early stages. At later 
times channeling of power from low |fc| to high |fc| is 
due to the steepening of structures by nonlinear effects 
(l4| : turbulence reaches fc-scales comparable to the 
inverse Debye length. Both these effects at high |fc| are 
electrostatic which is seen by comparing with Fig[3Jd); 
the magnetic field contains negligable power in deviation 
from that of a power-law on the inertial range. 

In summary, we have examined - numerically - tur- 
bulent wave spectra driven by the mixed mode instabil- 
ity (MMI) in 2D3V PIC simulations. Parameters for the 
simulations were chosen so as to favor the MMI over both 
the FI and the TSI. 

We interpret our main result as a consequence of weak 
turbulence, which transfers power turbulently on shells of 
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FIG. 3: (Color online) Power spectra for the total, P(\E x (k) +i E y (k)\ 2 ), the transverse, P(E x (k)), and the longitudinal, 
P(E y (k)), electric field components - (a), (b) and (c) respectively. The spectra are integrated azimuthally, (f> £ [0,7r/2], over 
shells in fc-space, {k x ,k y } = {|fc| cos <f), \k\ sin <j>} - cf. figlHb,c). Each curve correponds to one time instance, and these are 
distributed at equal intervals, t = [0, 10.., 100, 110]. Curves of low power (blue) gives the spectra for early times, and curves 
of high power (orange and yellow) reveils late time spectra, (d) Power spectrum of the total (out-of-plane) magnetic field, 
P(B z (k)), shown only for simulation end time, t = t en d = HO. The total magnetic-to-electric energy content at t en d is 
Eb/Ee ~ 10; a considerable amount of electromagnetic power is present but is still less than the electrostatic power. 



constant wavenumbers [lfj as seen in Fig^a-d). Neither 
the coalescence of filaments - transferring power to lower 
kr \l4 . 1 X 51 ] , nor phase space hole coalescence - transfer- 



ring power to lower k y 2, 17, 18(, however, can be cap- 
tured in an ideal weak turbulence model. This may ex- 
plain why we do not observe a neither simple ring distri- 
bution, nor constancy of the structure in Fig[TJa), along 
4> = const. 

Future work should be conducted to determine the rel- 
ative strength of the (ideally) weak and strong turbulence 
contributions to the total dynamics. In particular aspects 
of strong turbulence deserve thorough investigation. De- 
tailed analysis of the powerlaw slopes in FigJ3£a-d) is also 
important, for determination of the relative strengths 
of magneto-bremsstrahlung (synchrotron/jitter Q) and 
bremsstrahlung spectral contributions. 

A thorough parametric scan and a characterization of 
electromagnetic turbulence in shocks is important, when 
linking radiation spectra with predictions of physical 
conditions in a broad class of astrophysical jets. 
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